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A SET OF FORTRAN IV SUBROUTINES FOR GENERATING PRINTED PLOTS 

by Lois T. Dellner 
Lewis Research Center 

SUMMARY 

A set of subroutines, written in FORTRAN IV, easy for the FORTRAN programmer 
to use, provides printed plots as part of normal output. These subroutines are simpli- 
fied so that no choices must be made and generalized so that choices may be made by the 
programmer of the plotting characters, the scales, the appearance of the grid, and other 
options. The FORTRAN routines generate ordinary output records (up to 132 BCD char- 
acters) suitable for on- or off-line printing. 


INTRODUCTION 

Various versions (refs. 1 to 3) of this plotting system have been in use at Lewis 
Research Center since April 1962. The current version is a set of FORTRAN IV sub- 
routines which generates ordinary FORTRAN output records (up to 132 BCD characters) 
suitable for printing on- or off-line. 

The system permits, but does not require, the programmer to choose the plotting 
characters, the scales, the grid-line spacing, etc. It handles single or multiple curves, 
prints true scales, and permits titles above or below the plot. The next section is a 
manual for users. 

The routines are presently in use at Lewis on a 7094 Model II - 7044 direct-coupled 
system. All output is under the control of the 7044 supervisory system. The programs, 
however, are almost entirely machine independent, and the information provided in the 
section SYSTEM MANUAL is intended to make it very simple for the systems program- 
mer to make the changes required to adapt the plotting system to other machine configu- 
rations. The final section contains examples of plots obtainable with these subroutines 
(see pp. 61 to 65). 


USERS MANUAL 

This system offers printed plots as part of normal output with a minimum of pro- 
gramming effort. The programmer writes CALL PIXDTXY (for a single curve) or CALL 
PLOTMY (for multiple curves). The arguments, or call list, include the names of the 
arrays to be plotted and specify the number of points per curve and the number of curves. 
He may precede the call by writing a title to be printed above the plot. He may follow the 
call by writing a legend to be printed at the bottom of the plot. 

The plot or plots are printed as part of the regular output listing with no delay. No 
changes in his card-handling procedures nor special instructions for the operators are 
required. 

If he is using PLOTXY, the values of the variable to be plotted in the x-direction 
must be in sequence. If they are not, the subroutine SORTXY (which makes the neces- 
sary rearrangements) is supplied to be used before calling PLOTXY. For either PLOTXY 
or PLOTMY, if the size of the elements in (or the total range of) any array is not known 
to be within certain limits, the programmer calls the subroutine SCALE for each array 
before calling the plotting subroutine. SCALE will transform the array to suit PLOTXY 
and PLOTMY only if it is necessary. 

In addition to the minimum -effort use just described, the programmer may choose to 
use one or more of several options that permit him to control, for example, the appear- 
ance of the grid (by specifying the frequency of the grid-lines in either direction), the 
scale for either variable (by specifying the scale factor and a starting value), the plotting 
character, etc. 

This section contains detailed instructions for the use of PLOTXY and PLOTMY and 
brief descriptions of SCALE and SORTXY. Some assistance in debugging users pro- 
grams completes this section. 

I. PLOTXY 

To get plotted output using PLOTXY, the corresponding pairs of ordinates to be 
plotted must be in two arrays. For an example, let us name the arrays XDOWN and 
YACROS and assume each is NPTS elements long. These names are chosen specifically 
to call the user's attention to the fact that the x-direction is down the page. 

A. Call 

Write CALL PLOTXY (XDOWN, YACROS, KODE, P) 

XDOWN is the name of the array containing the values of the variable to be plotted 
on the X scale (down the page) . The elements in the array are restricted 
as follows: 
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(1) They must be in floating point. 


(2) The absolute value of each element must^ be within permissible 
limits (approximately 10" < e < 10°). 

2 

(3) They must be in order, either increasing or decreasing. 

YACROS is the name of the variable to be plotted on the y scale (across the page) . 
The elements of the array are restricted exactly as in (1) and (2) above. 

KODE is the name of an integer . Many options are provided for the user of 

PLOTXY. Each has a number associated with it. The sum of the num- 
bers representing the options being used is KODE (0 < KODE < 127) . For 
your first plot, use KODE = 0. When KODE = 0, the starting-values and 
scale-factors in both directions are computed by PIOTXY. The other 
effects of KODE = 0 are shown in the example plot; 

(1) The plotting character is an asterisk. 

(2) The grid-line frequency is 10 x 10. 

(3) Nothing is printed to the left of the plot. 

P is an array. 

P(l) must contain NPTS (the number of points to be plotted) in floating 
point. 

If KODE = 0 or 64, there are no other requirements for P. 

For all other values of KODE, the requirements for the P array are dis- 
played (see section D. Using the Options, p. 6). 


^If the element size of any array is unknown or out of range, write CALL SCALE 
(NPTS, A, KA) before calling PLOTXY (see HI, p. 13). 

^If the array to be plotted on the x scale is not in order, write CALL SORTXY 
(XDOWN, YACROS, NPTS) before calling PIX5TXY (see IV, p. 13). 
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B. Title 


The call for PIX)TXY may be preceded by writing at least one TITLE (PT) record. 
A representative example is: 

WRITE (6, 500) 

500 FORMAT (2HPT, 71X, 14HSAMPLE EXAMPLE/2HPT, 69X, 

17H JOLO SYSTEM PLOT) 

The result is shown on the plot on the facing page. 

C. Legend 

The call for PLOTXY may be followed by writing at least one LEGEND (PL) 
record. A representative example is: 

WRITE (6, 502) 

502 FORMAT (2HPL, 70X, 16HF SUB2 VS. THETA/2HPL, 70X, 

16H(— ) (RAD)) 

The result is shown on the plot on the facing page. 
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SAMPLE EXAMPLE 
JULO SYSTEM PL3T 
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SUB2 VS. THETA 



D. Using the Options 

Each option has a number associated with it. The sum of the numbers repre- 
senting options chosen is KODE, the third argument in the call. 


TO CHOOSE: \ 


AND 

SUPPLY 

IN 

(FORMAT) 

KODE 0 
USES 

The plotting character. 

1 

Desired 

P(2) 

1H0 

♦ 

0 represents any acceptable FORTRAN 


charac- 




character except the minus sign. 


ter 




The frequency of x grid>lines. 

2 

m 

P(3) 

Floating 

10 

They are printed every m line-spaces 
below the first. If m = 0, only two 
will be printed, one above and one 
below the plot. 




point 


The frequency of y grid-lines. 

4 

n 

P(4) 

Floating 

10 

They are printed every n positions 
to the right of the first. If n = 0, 
only one is printed, at the left of 
the plot. 




point 


The X scale. ^ 

16 

Three 

P(6), 

Floating 

Scale 

The scaling parameters FX, repre- 


scaling 

P(7), 

point 

computed 

senting the starting- value , and DX, 


param- 

P(8) 


by 

representing the scale -factor for one 


eters, 



plotting 

line space, must be whole numbers of 


KSX, FX, 



routine 

magnitude less than 10®. (0 < KSX < 6) 
Seel.-E. Scaling. 


andDX 




The y scale. ^ 






The scaling parameters FY, repre- 

32 

Three 

P(9). 

Floating 

Scale 

senting the starting-value, and DY, 


scaling 

P(10), 

point 

computed 

representing the scale -factor for one 


param- 

P(ll) 


by 

print position, must be whole numbers 


eters. 



plotting 

of magnitude less than 10®. (0 < KSY < 6) 


KSY, FY, 



routine 

See I. -E. Scaling. 


andDY 




To print P array. 

8 

Array 

pa + 11) 

BCD 

No 

The 6 BCD characters in P(I + 11) will 


to be 

[through 


print- 

be printed at the left of the plot on the 
line on which the point (x.,y.) is plotted. 
(If desired material is numeric, write it 
out using format F6. n, and read it back 
into P array using A6) 


printed j 

lp(NPTS + 11) 


out 

To print coordinates. 

64 

No re- 



No 

The coordinates of each point will be 


quire- 



print- 

printed on the line on which the point is 
plotted. 


ments 



out 


*To specify only FX, omit DX in P(8) and place 1. in P(5). 

To specify only FY, omit DY in P(ll) and place 2. in P(5). 

To specify FX and FY only , omit DX and DY and place 3. in P(5). 
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E. Scaling 


When the user wishes to specify his own scale in either direction, in addition to 
increasing KODE by 16 or 32, he must: 

(1) Choose his desired starting-value F. 

(2) Choose his desired scale-factor D - for one line-space if he is specifying the 

X scale, for one print position if he is specifying the y scale. 

(3) Determine a value of N such that: 

N 

(a) F X 10 is a whole number. 

N 

(b) D X 10 is a whole number. 

(c) 0<N<6. 

(4) Compute 6 - N. 

The integer calculated in step (4) is the first of three scaling parameters that must 
be available in the P array (in P(6) if scaling x, P(9) if scaling y) when the plotting sub- 
routine is called. F x 10^ is the second scaling parameter, and is placed in P(7) or 
P(10). F may be zero. D x 10^ is the third, and is placed in P(8) or P(ll). Note that 
D may never be zero, and when scaling x for PLOTMY, D must be positive. 

If the user wishes to supply only the starting value and leave the choice of scale- 
factor to the plotting system, he must do steps (1), (3) (a), (3)(c), and (4) above, and 

(a) Supply 1. in P(5) if option 16. 

(b) Supply 2. in P(5) if option 32. 

(c) Supply 3. in P(5) if both 16 and 32 are being used. 
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n. PLOTMY 


Although more than one curve can be plotted on the same grid with PLOTXY, only 
one plotting character will be used for all curves. PIX)TMY provides a different plotting 
character for each curve. Options, similar to those in PLOTXY, are available but not 
required, except that option 1 must be used if there are more than six curves. However, 
a selection must be made from three Variations: DUPX - when more than one set of y 
values corresponds to the same set of x values; DUPY - when more than one set of x 
values corresponds to the same set of y values (this offers complete control of which 
variable is to be plotted in which direction) ; and NO DUP - when each set of x values has 
a corresponding set of y values. 

A. Call 

Write CALL PLOTMY (XDOWN, YACROS, KKK, P) 

XDOWN is the name of the array containing the values of the variable to be 

plotted on the x scale (down the page) . The minimum DIMENSION of 
this array depends on the Variation selected. (See E. Variation Layout, 
pp. 10 to 12). The elements of this array are destroyed by PLOTMY. 
The elements are restricted as follows: 

(1) They must be in floating point. 

(2) They must be within permissible range. (See I. PLOTXY, p. 3.) 

YACROS is the name of the array containing the values of the variable to be plot- 
ted on the y scale (across the page). The minimum DIMENSION of this 
array depends on the Variation selected. (See E . Variation Layout. ) 

The elements of this array are restricted exactly as in (1) and (2) above. 

KKK is the name of an array. The first element must be KODE (as in 

PLOTXY, this is the sum of the option numbers), and the second ele- 
ment must be KN (the number of curves). The third element must be: 

(1) The number of points in one curve for DUPX and DUPY. 

(2) The number of points in the first curve for NO DUP. 

The remaining odd-numbered elements are only required for NO DUP, 
and the remaining even-numbered elements are only required for op- 
tion 1. (See E. Variation Layout.) The minimum DIMENSION of KKK 
must be 2 * KN + 2 or 14, whichever is greater. 
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P is the name of an array. The first element must specify the Variation 

selected: 


P(l) = 1. for DUPX 
P(l) = 3. for DUPY 
P(l) = 5. for NO DUP 

The remaining elements and the minimum DIMENSION of P are func- 
tions only of the options being used. (See E. Variation Layout.) 

B. Title 

A TITLE record may be written before calling PIXJTMY, exactly as for PIX)TXY 
(p. 4). 

C. Legend 

A LEGEND record may be written after calling PLOTMY, exactly as for PIOTXY 
(p. 4). 

D. Using the Options 

The use of the options is exactly the same as for PLOTXY (see I. -D. Using the 
Options, p. 6) with the following exceptions: 

(1) If option 16 is used, DX must be positive. 

(2) When option 64 is used, only the ordinate in the x-direction will be printed. 

(3) Unless option 1 is chosen, the plotting characters *, +, 0, X, =, 0 will be 
used for the first six curves. To use others, use option 1 and supply the desired 
plotting character for the first curve in KKK(4), for the second curve in 
KKK(6), .... If more than six curves are to be plotted, option 1 must be used, 
and all plotting characters supplied. To put a plotting character, for example, $, 
in KKK(4) , write 

DATA KKK(4)/1H$/ 
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KKK 


KODE 

KN 

NPTS 

[PCl]j 

[PC2]^ 

[PC3]i 

[PC4]j 


See 


below 






DUPY 


[ ]j^ means supply item in brackets if option N is being used. 



DIMENSIONS 


X(KN * NPTS) 

Y(KN ♦ NPTS) 

KKK = (2 * KN + 2) or (14), whichever is greater. 

If option 8 is not used, P(<11), depending on 
use of options 2, 4, 16, and 32. 

If option 8 is used. 


P(KN * NPTS + 11) 


NO DUP 


[ ] means supply item in brackets if option N is being used. 



11 + NPl + NP2 
+ NP3 + . . . 


DIMENSIONS 

X(NP1 + NP2 +NP3 + . , 

Y(NP1 +NP2 +NP3 + . . .) 

KKK = (2 ♦ KN + 2) or (14), whichever is greater 

If option 8 is not used, P(<11), depending on 
use of options 2, 4, 16, and 32. 

If option 8 is used, 

P(ll + NPl + NP2 + NP3 + . . .) 


m. SCALE 


CALL SCALE (NPTS, A, KRSTR) 

This subroutine finds the largest absolute value of the NPTS elements of A and com- 
putes the characteristic of its logjQ using a function subprogram KHAR (described under 
PISTUG, p. 46). If the characteristic K is -2 <K < 4, KRSTR is set to zero and control 
returns to the calling program. If K > 4 or < -2, each element of A is multiplied by a 
power (KRSTR) of 10 to transform the array to suit PLOTXY and PLOTMY. KRSTR is 
returned to enable the user to ReSToRe the array or to record how it has been altered. 

IV. SORTXY 


CALL SORTXY (V, W, NPTS) 

This subroutine rearranges the NPTS elements of the V array in order of increasing 
size. The elements of the W array are moved to maintain the original pair -relation; that 
is, if the fifth element of the V array is moved to the first position of V, the fifth element 
of W is moved to the first position of W. 

V. ERROR MESSAGES 

A. The message BAD LABELS printed below and to the left of the plot indicates that 
one or more of the x or y grid-labels may be incorrect. This can be caused by DX 
or DY too large or too small or by a label requiring more than the allotted nine print 
positions. 

B. When the system computes the scale factor for either array and finds the range is 
zero, calculations are stopped. For either PLOTXY or PLOTMY the plot terminates 
with a message : 

(1) For PLOTXY - N.G. followed by the contents of X(l), Y(l), X(2), Y(2), KODE, 
and P(l) 


or 

(2) For PLOTMY - N.G. followed by the contents of X(l), Y(l), X(2), Y(2), K(l), 
K(2), andK(3). 


C. If the values in the X array are not in order when PLOTXY is called, or if you use 
option 16 and a value of X lies outside your specified starting value, the message 

X OUT OF ORDER. I=NNNNN 

will be printed followed by the message described in (2) on the previous page, and the 
plot is terminated. 

D. In PLOTMY if the plotting character search fails, the message ERROR in K ARRAY 
is printed. This usually means DIMENSION of KKK is incorrect. 


SYSTEM MANUAL 

The objective of this section is to simplify the work of the systems programmer who 
implements and may have to modify this plotting system. To such a programmer the sys- 
tem consists of a set of FORTRAN IV subroutines (PLOTXY, PLOTMY, PISTUG, 

SORTXY, SCALE) which are executed in a machine configuration that controls all input- 
output. A brief description and a complete listing are included for all routines. For 
PLOTXY, PLOTMY, and PISTUG, comprehensive block diagrams and dictionaries of the 
FORTRAN variable names are supplied. These are intended as systems debugging aids . 
Each dictionary is designed for use primarily as a cross-reference with the corresponding 
FORTRAN listing and block diagram. Within each dictionary, each variable name is iden- 
tified in terms of its contribution to the program. Throughout this section, the subrou- 
tine name is printed at the bottom of each page for assistance in rapid cross-referencing. 

The final discussion in this section is of the input-output structure within which the 
system operates at Lewis. 

I. PLOTXY 

A. Description 

To plot one curve, the programmer will usually use PLOTXY. His call statement 
lists the names of the arrays to be plotted, the number KODE (an indicator of which 
options are being used), and finally, an array whose first element is the number of 
points to be plotted and whose other elements are any additional data required by the 
choice of option. 

PLOTXY writes a blank line (FORMAT 2HPT; see VI, p. 58) to initiate the plot, 
and calls the auxiliary subrouting PISTUG (p. 46) to compute scaling parameters for 
either or both scales and check label sizes. Communication with PISTUG is through 

COMMON/JOLO/N, F, DX, XYX, FORY, STUG, LABOUT, TONLY,KSW64,KPWR,KFD, TLINX 
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PLOTXY 



Records are then written for the plot beginning with a line of 11 y grid- labels. 

For each line of the plot any required horizontal and vertical grid-line characters as 
well as correctly positioned plotting characters for all points in that line are written. 

An X grid-label is printed every tenth line using space suppress. When all points have 
been plotted, the plot is terminated (with an x grid-line at the next x grid-label, fol- 
lowed by a line of y-labels and a final blank line (FORMAT 2HPL; see VI, p. 58)). 

Since the values in the X array are required to be in order (up or down), no sorting 
is required and the user's arrays are undisturbed. 


PLOTXY 
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B. Program Listing 


SUBROUTINE PLOTXY(X,Y»K,P) 

COMMON/JOLO/N.F»OX,XYX,FORY,STUG,LABOUT,TONLY,KSW64,KPWR.KFO,TLINX 
LOGICAL XYX.FORY. STUG»TONLY,XGL,LS 
DIMENSION xm«Y{ll«P(l) 

DIMENSION FLSI3).FLABI4) ,FYLA8(4) .YLABELdl ) tA(l04) tELSO) 
EQUIVALENCE ( FLAB ( 3 1 , IFLAB3) , C FYLAB (3 I » I YLAB) 

DATA MASK!. MASK2,MASK4*MASK8»MASK16tMASK32 tMASK64 / 

1 01* 02* 04* 010* 020* 040* 0100 / 

DATA FYLABf 1)*FYLABI2I*FYLA83*FYLAB(4) / 

I 6HC2HP **6H20X*11 *6H F10.0*1H) / 

DATA BLANK*XGRID*YGRID /IH ,1H-,1H1 / 

DATA PCSTD/IH*/ 

DATA RMARK/0726060606060/ 

DATA FLS( 1)*FLS28*FLS264*FLS2B*FLS38*FLS3B / 

1 6HI 2HP+«*4H12X* * 3H6X * *3HA6* *3HA6) *6H2F6.3 I / 

DATA FLABI 1)*FLAB(2)*FLAB3*FLAB(4) / 

1 6H(2HP4-**4H18X**6H F9.0*1H) / 

100 WRITE (6*5001 
500 FORMAT! 2HPT1 
102 KODE=K 
N*P( 1) 

LABOUT*! 

FLABI 3) * FLAB3 
LS = .FALSE. 

FYLAB( 3)=FYLAB3 
KSWBsO 
KSW64=0 
110 PC=PCSTO 

112 IFUAND(KOOE*MASKII i.GT.O. }PC « P(2) 

114 M=10 

116 IF( ( AN0(K0DE*MASK2n .GT.O. ) M-P( 3) 

117 IF (H.EO.OIM^IOOO 

118 NY*10 

120 IF( I AN0(K0DE*MASK4H .GT.O. »NY = PI4) 

121 IFCNY.EO.O) NY*1000 

122 IF( ( AN0(K00E*HASK64n.GT.0.)KSW64-2 

124 IF( (AND(K00E*MASK8) I.GT.O. 1KSW8-1 

125 K864>KSW8+KSW64 

126 lF(K864-2) 132* 128*138 
128 FLS( 2)aFLS264 

130 60 TO 139 
132 FLS(2»*FLS28 
134 FLS(3)=FLS38 
136 60 TO 140 

138 FLS(2»=FLS2B 

139 FLS(3)*FLS38 
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140 XYX*. FALSE. 

142 F0RY=.TRUE. 

144 STUG=. FALSE. 

146 TONLY*. FALSE. 

C 

148 IFC (AN0(K00E«MASK32).LE.0.n GO TO 172 

151 STUG=.TRUE. 

152 KSY=P(9> 

154 PWRl0Y = 10.**(KSY-6) 

156 FY »P( 10)*PHR10Y 
158 F * FY 
C 

160 IF(P( 5I.GE.2.) GO TO 172 
162 TONLY=.TRUE. 

164 0Y= P< 11 )»PWR10Y 
166 DX= DY 

172 CALL PISTUG(Y) 

173 IFIDX.EO.O.) GO TO 700 

174 FY*F 
176 DY=DX 

180 IF(KSW64.E0.2) KPWRY»KPUR 
190 IYLAB=IYLAB-KFO 
C 

200 XYX *.TRUE. 

202 FOR Y=. FALSE. 

204 STUG*. FALSE. 

206 TONLY=. FALSE. 

208 TLIMX=55*( l+N/35) 

C 

210 1FC(ANO(KOOE,MASK16).LE.O. n GO TO 232 

213 STUG=.TRUE. 

214 KSX - P(6) 

216 PMR10X*010.**IKSX-6J 
218 FX= PI 7)4PWR10X 
220 F=FX 
C 

222 IFIHODI IFIXIPI 5)).2) .EO.ll GO TO 232 
224 TONLY».TRUE. 

226 OX =P(8)«PWR10X 
232 CALL PISTUGIXI 

IFIDX.EO.O. I GO TO 700 
234 FX*F 

240 IFIKSW64.E0.2I KPWRX=KPWR 
248 IFLAB3*IFLAB3-KF0 
C 

250 IFIKSW64.EO.O)GO TO 264 

252 KOUTX=-KPWRX 

254 KOUTY=-KPWRY 

256 F10X*10.**KPWRX 

258 F10Y=10.**KPWRY 

260 WRITE 16.502) KOUTX.KOUTY 

502 FORMAT! 2HPT.7X. 3HX4E,I2,4H Y*E.I2) 

C 

264 DO 278 1=1,11 

266 TEMP = FY+FL0ATII-l)*DY*10. 

268 ATEMP= ABSITEMP) 


PLOTXY 
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270 IF < ATEMP.lt. 1.E-7I TEMP * 0. 

272 IF f ATEMP.6E.l.E>7)LABOUT»2 
278 YLABELI I»=TEMP 
300 KSYLAB =l 

302 WRITE (6.FYLABI I YLABEL 41 I ,I ,11 ) 
304 60 TO 4 306, 7001, KSYLAB 
C 

306 KSYLAB *2 
310 LCTR*0 
NCTR»1 
K0UT»1 
KOUIT* 1 
C 

320 IF4M004LCTR*Mn328,322,328 
322 AFILL» XGRIO 
324 XGL ».TRUE. 

GO TO 330 
328 XGL =. FALSE. 

AFILL»BLANK 
330 DO 332 I«2,104 
332 A4 I) s AFILL 
334 DO 336 I»2,104,NY 

336 A4 I )- Y6RI0 
Al 1) sBLANK 

338 GO TO 4 340.4001 ,KOUT 
C 

340 KX »4 X4NCTR )-FXI/DX +.5 
342 IF4KX-LCTRI630.350.600 

350 KY* 4 Y4NCTR l“FY»/0Y+.5 

351 LS* .TRUE. 

352 TPC * PC 

353 KYL * KY+2 

354 IF4KY.lt. 0) GO TO 360 
356 IP4KY.GT.101)G0 TO 364 
358 GO TO 370 

360 KYL*1 
362 60 TO 366 
364 KYL»104 
366 TPC»RMARK 
C 

370 A4KYLI «TPC 
372 J*1 

374 IF4KSW8.EO.O) GO TO 380 
376 ELS4JI»P4NCTR+11) 

378 J= J+1 

380 IF4KSW64.E0.0) GO TO 386 
382 ELS4J|sX4NCTR)/F10X 
384 ELS4J + 1I * Y4NCTRI/F10Y 
C 

386 IF4NCTR.6E.NIGO TO 392 
388 NCTR4NCTR4-1 
390 GO TO 340 
C 

392 KOUT * 2 
394 M* 10 
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400 IF! XGL .AND. fMQOtLCTR , 10) .EO.OI )KOUI T»2 
C 

600 WRITE (6.504)IAf n,I»l,104) 

504 FORMAT! 2HP «26X*104A1) 

602 IFfK864.NE^0) GO TO 620 
C 

604 IF(HOOILCTR.IO) )614«606t614 
606 XLABEL *FX+FLOAT(LCTR )*0X 
608 TEMP »ABSI XLABEL) 

610 IFCTEMP.GE.1.E+7)LAB0UT*2 

612 IFITEMP.lt. l.E-7)XLABEL=0. 

613 WRITE! 6. FLAB) XLABEL 
C 

614 LCTR » LCTR+1 

616 GO TO ! 320,302)*K0UIT 
C 

620 IF!.NOT.LS) GO TO 604 
622 LS s .FALSE. 

624 WRITE !6.FLS)!ELS!I ).I»1.K864) 

626 GO TO 604 
630 L ABOUT* 4 
C 

700 GO TO ! 710,702.706,704),LABOUT 
702 WRITE !6, 506) 

506 F0RMAT!2HPL.3X«10HBA0 LABELS 
GO TO 720 
C 

704 WRITE! 6. 520) NCTR 

520 FORMAT! 2HPL.18HX OUT OF ORDER. I».I5) 
706 WRITE!6. 508)!X! I ) * Y! I ) *I *1 *2) « K* P! 1 ) 
508 FORMAT! 2HPL.5H N.6. «4G20. 8«I 6«f8.2) 

708 GO TO 720 
710 WRITE! 6. 510) 

510 FORMAT !2HPL) 

720 RETURN 
END 
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C. Block Diagram 
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D. PLOTXY Dictionary 


AFILL 

ATEMP 

BLANK 

DX 

DY 

ELS 

F 

FlOX 

FlOY 


Array in which each line of the plot proper is constructed. It is ini- 
tialized to 104 blanks or 104 minuses (the horizontal grid-line charac- 
ter). The vertical grid-line characters (I's) are placed in position 2 
and every NY**^ position thereafter. Finally, the plotting characters 
(if any) are positioned. 

Temporary storage for the character (1 or -) with which the A array 
is initialized. 

Temporary storage for the absolute value of each y grid-label. 

BCD representation for a blank printing character. It is used to ini- 
tialize the A array for lines that are not x grid-lines. 

Vertical scale factor (per printing line) . It is either calculated by 
PISTUG based on the range and length of the X array or it is calculated 
from values supplied by the user of option 16 in P(6) and P(8) as 
DX = P(8)- Also see DX, COMMON/JOLO/ - Dictionary, 
page 53. 

Horizontal scale factor (per print position) . It is either calculated by 
PISTUG based on the range and length of the Y array, or it is calcu- 
lated from values supplied by the user of option 32 in P(9) and P(ll) as 
DY = P(ll)-10^^®^“®. Also see DX, COMMON/JOLO/ - Dictionary, 
page 53. 

Array in which the information (the first 18 characters of each line) to 
be printed at the left of the plot is assembled. See FLS, page 24. 

See COMMON/JOLO/ - Dictionary, page 53. 

IqKPWRX divisor of all values in the X array before they 

are printed at the left of the plot when option 64 is used. 

IqKPWRY jg divisor of all values in the Y array before they 

are printed at the left of the plot when option 64 is used. 
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FLAB 

Variable FORMAT array with which each x grid-label is written. It 
is of the form 

(2HP+, 18X, F9.d) 

where 0 < d < 6 . A value of d is inserted by the program (as KFD) 
after the call to PISTUG to process the X array. See KFD, page 53. 

FLAB3 

Initial value for FLAB (3) . 

FLS 

Variable FORMAT array with which the output to the left of the x 
grid-labels is written out if option 8 or option 64 is used. It is con- 
structed during execution and becomes one of three different state- 

Used to write: 

If 8 only (2HP+, 12X, A6) One word of P array 

If 64 only (2HP+, 6X, 2F6.3) The modified coordinates 

If 8 and 64 (2HP+, A6, 2F6. 3) One word of P array and the 

modified coordinates 

FLS264 

Used to fill FI£(2) if option 64. 

FLS28 

Used to fill FLS(2) if option 8. 

FLS2B 

Used to fill FLS(2) if both option 8 and option 64 are used. 

FI538 

Used to fill FLS(3) if option 8. 

FLS3B 

Used to fill FLS(3) if option 64 or both 8 and 64 are used. 

FORY 

See COMMON/JOLO/ - Dictionary, page 53. 

FX 

Starting value for the vertical scale. It is either calculated by PISTUG 
based on the range, length, and XI (see p. 52) of the X array, or it is 
calculated from values supplied by the user of option 16 in P(6) and 
P(7) as FX = P(7)- 

FY 

Starting value for the horizontal scale. It is either calculated by 
PISTUG based on the range, length, and XI (see p. 52) of the Y array, 
or it is calculated from values supplied by the user of option 32 in 
P(9) and P(10) as FY = P(10)- 
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FYLAB Variable FORMAT array with which the line of y grid- labels is 
written. It is of the form 

(2HP , 20X, llFlO.d) 

where 0 < d < 6. A value of d is inserted by the program (as KFD) 
after the call PISTUG to process the Y array. See KFD, page 53. 

FYLAB3 Initial value for FYLAB(3) . 

I Index of loop forming y grid-labels. 

IFLAB3 Equivalent to FLAB(3). 

lYLAB Equivalent to FYLAB(3) . 

J Subscript used when filling ELS array. 

K864 Sum of KSW8 and KSW64. It has the possible values 0, 1, 2, 3. 

KFD See COMMON/JOLO/ - Dictionary, page 53. 

KHAR Function subprogram described under PISTUG, page 46. 

KODE Sum of the option numbers. It is supplied by the programmer in the 
third argument of the call. Branching on single bit positions using 
MASKl, MASK2, MASK4, MASK8, MASK16, MASK32, and MASK64 
serves to identify the options being used. 

KOUT Switch initialized to 1 and set to 2 after all points have been processed. 

KOUTX Value in the special heading of option 64 for the x coordinates. It is 
the negative of KPWRX. 

KOUTY Value in the special heading of option 64 for the y coordinates. It is 
the negative of KPWRY. 

KPWR See COMMON/JOLO/ - Dictionary, page 53. 

KPWRX KPWR for the X array. 

KPWRY KPWR for the Y array. 

KQUIT Switch initialized to 1. It is set to 2 just before the last line (a labeled 
X grid-line) of the plot is printed. 
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KSW8 

KSW64 

KSX 

KSY 

KSYLAB 

KX 

KY 

KYL 

LABOUT 

LCTR 

LS 

M 


Switch initialized to zero and set to 1 if option 8 is used. 

See COMMON/JOLO/ - Dictionary, page 53. 

Fixed point form of the value supplied in P(6) by the programmer using 
option 16. P(7) and P(8) are multiplied by to get FX and 

DX. 

Fixed point form of the value supplied in P(9) by the programmer using 
option 32. P(10) and P(ll) are multiplied by to get FY and 

DY. 


Switch initialized to 1 and set to 2 after the first line of y grid-labels 
has been written. 


Line number of a particular point. It is computed from 
rounded. The first x grid-line is line number zero. 


[(x. - FX)/DX] 


’’Printing position” of a particular point. It is computed from 
(y^ - FY)/DT2^ rounded. The first y grid-line is in ’’printing posi- 
;ion” zero - actually the second element of the A array. 

Value of KY + 2 if the y-coordinate falls on the plot; 1 if KY is 
negative; 104 if KY > 101. 

See COMMON/JOLO/ - Dictionary, page 53. 

Line counter. It is set to zero for the first x grid-line and is 
stepped by 1 after each line is printed. 

Logical variable set to TRUE whenever a point has been found for the 
current line. If option 8 or 64 are in effect, LS is used to control a 
branch (to print on the left side of the plot) and is set to FALSE imme- 
diately after the branch. 

Frequency of the x grid-lines. It is either supplied in P(3) by the 
programmer if option 2 is used or set equal to 10. If the programmer 
specifies zero, M is set to 1000 and only the first and last x grid- 
lines will be printed. 
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MASKl 

MASK2 

MASK4 

MASKS 

MASK16 

MASK32 

MASK64 

N 

NCTR 

NY 


P 

PC 

PCSTD 

PWRIOX 

PWRIOY 

RMARK 

STUG 

TEMP 




VSee KODE, page 25. 


See COMMON/JOLO/ - Dictionary, page 54. 

Counter for the number of points that have been processed. It is ini- 
tialized to 1 and tested for equality to N. It is increased by 1 each 
time the test fails. 

Frequency of the y grid-lines. It is either set to 10 or supplied in 
P(4) by the programmer using option 4, If it is supplied as zero, 
it is replaced by 1000 and only the left-most y grid-line will be 
printed. 

Array name, the fourth argument in the call list of this subroutine. 
P(l) contains number of points to be plotted. Further contents are 
prescribed by the options being used. 

Plotting character. It is either supplied by the programmer using 
option 1 in P(2) or by the program as PCSTD. 

Standard plotting character * . 

Value of 10*^^"®. 

Value of 10*^®^"®. 

Record mark character (4^). 

See COMMON/JOLO/ - Dictionary, page 54. 

Temporary storage. 
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TUNX 

TONLY 

TPC 

X 

XGL 

XGRID 

XLABEL 

XYX 

Y 

YGRID 

YLABEL 


See COMMON/JOLO/ - Dictionary, page 55. 

See COMMON/JOLO/ - Dictionary, page 55. 

Current plotting character. 

Array name, the first argument in the call list of this subroutine, 
containing the values of the variable to be plotted on the vertical scale 
(down the page). They must be in sequence either increasing or de- 
creasing. 

Logical variable set to TRUE whenever the value of LCTR is evenly 
divisible by M. 

X grid-line character, a minus sign. 

Current value of the x label. It is computed only for every tenth line 
from FX + LCTR-DX. 

See COMMON/JOLO/ - Dictionary, page 55. 

Array name, the second argument in the call list of this subroutine, 
containing the values of the variable to be plotted on the horizontal 
scale (across the page) . 

y grid- line character - the digit 1. 

Array holding the 11 labels for the y axis. 


28 


PLOTXY 


n. PLOTMY 


A. Description 

Although more than one curve can be plotted on the same grid using PLOTXY, they 
will all use the same plotting character. When a different plotting character for each 
curve is desired, the programmer uses PIX)TMY. The CALL statement lists the 
names of the arrays to be plotted; an array containing the number KODE, the number 
of curves, and possibly other information - depending on the options used and the 
Variation chosen; and finally an array whose first element indicates the Variation 
chosen and whose other elements are any additional data required for the options being 
used. 

The records written by PLOTMY are identical to those written by PLOTXY 
(p. 14), and the auxiliary subroutine PISTUG is used in the same way. 

Because the X array is not required to be in order, a search for the next point to 
be plotted occurs for every point. The search technique used avoids the use of addi- 
tional storage by using the X array itself in such a fashion that the values in the array 
are destroyed. The Y, KKK, and P arrays are undisturbed. 
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B. Program Listing 


SUBROUTINE PLOTMVIX,Y.K,P» 

C 

COMMON/JOLO/N,F*OX.XyXfFORY»STUGtLABOUTtTONLY»KSW64«KPWRfKFOfTLINX 
LOGICAL XYX.FORY, STUG.TONLY, XGL «LS 
DIMENSION xm,Y(l)tP(l)fKll) 

DIMENSION FLS(3)«FLAB(4) »FYLAB<4) * YLABEL ( 11 ) «A( 104) ,ELS(3) 
DIMENSION KPCST0I6) 

EQUIVALENCE ( FLAB ( 3 ) • IFLAB3) . (F YLAB 13 ) ,1 YLAB) 

EQUIVALENCE IKPC.TPC) 

DATA MASKl, MASK2*MASK4«MASK8*MASK16 «HASK32 .MASK64 / 

1 01* 02* 04* 010* 020* 040* 0100 / 

DATA FYLABI1)*FYLAB(2)*FYLAB3*FYLAB(4I / 

1 6HC2HP **6H20X*11 *6H F10.0*1H) / 

DATA BLANK*XGRID*YGRID /IH *1H-*1H1 / 

DATA RMARK/072 6060606060/ 

OATAIKPCSTDI I) *I»1*6)/1H«*1H4^ *lHO*lHX*lHs*lHO/ 

DATA FLS< 1)*FLS26*FLS264*FLS2B*FLS38*FLS3B / 

1 6H(2HP-f**4H12X**4H12X**4HA12**3HA6) *6H F6.3) / 

DATA FLABI 1)*FLABI2)*FLAB3*FLABI4) / 

1 6H( 2HP4-**4H18X**6H F9.0*1H) / 

C 

100 WRITE 16*500) 

102 KODE=K( 1 ) 

KN-KI 2) 

NPTS*KI 3) 

LAB0UT«1 
FLABI 3) - FLAB3 
LS = .FALSE. 

FYLABI 3)*FYLAB3 
KSWa»0 
KSW64»0 
KTL»1 

KSWIsl 

KSWII^l 

11000 IF(P(l)-2.5) 11002* 11002*11010 
11002 KSUIs2 
11004 KTL^MCN 
11006 GO TO 11034 

11010 IFIPI l)-4. 0)11020*11020*11012 

11012 NPTST=0 

11014 DO 11016 I»1*KN 

11016 NPTST*NPTST+KI2*I+1) 

11017 TLINXs55«Il+NPTST/f 354KN) ) 

11018 GO T0112 

11020 KSWII«2 

11021 KTIMES»KN-1 

11022 00 11032 I*1*KTIMES 
11024 MM=I*NPTS 

11026 KI 2*I + 3)*NPTS 
11028 DO 11032 II»1*NPTS 
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11030 L»MM^II 
11032 Y«L1*YIII) 

11034 NPTST=KN*NPTS 
11036 TLINX*55*Cl+NPTS/35> 

112 IFI (AN0(K00E,MASKin. GT. 0.160 TO 114 

11112 DO 11113 I»1,KN 

11113 K(24l+2»*KPCSTDn » 

114 M»10 

116 IFI (ANO(KOOE.MASK2)).GT.O. )H«P(3I 

117 IF (H.E0.0)M»1000 

118 NY»10 

120 IFI IAN0IK0DE.MASK4I I.GT.O. )NY » PI4) 

121 IFINY.EO.OI NY»1000 

122 IFI IAN0IK00E«MASK64n.6T.0.)KSW64s2 

124 IFI lANOIKOOE. MASKS)). GT. 0. )KSW8sl 
C 

125 K864-KSW8+KSW64 

126 IFIK864-2) 132*128, )38 
128 FLSI2)sFLS264 

130 GO TO 139 
132 FLS{2)«FLS28 
134 FLSI3)«FLS38 
136 GO TO 140 

138 FLSI2)»FLS2B 

139 FLSI3)»FLS3B 
C 

140 XYXs. FALSE. 

142 FORY*.TRUE. 

144 STUGs. FALSE. 

146 TONLY*. FALSE. 

C 

148 IFI IANOIKODE,HASK32).LE.O. ) ) GO TO 172 
151 STUG*.TRUE. 

15? KSY=PI 9) 

154 PWRIOY^IO.** IKSY-6) 

156 FY »PI 10)*PWR10Y 
158 F s FY 
C 

160 IFIPI5).GE.2.) 60 TO 172 
162 TONLYa.TRUE. 

164 OYa PI 11)*PWR10Y 
166 DXa DY 
C 

172 NaNPTST 
11172 CALL PISTU6IY) 

173 IFIOX.EO.O.) GO TO 700 

174 FY=F 
176 DYaOX 

190 IYLAB»IYLAB-KFO 
C 

200 XYX=. FALSE. 

202 FORY*. FALSE. 

204 STUG*. FALSE. 

206 TONLY*. FALSE. 

C 

210 IFI IAN0IK00E,MASK16).LE.0.) ) GO TO 232 
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213 STUG=.TRUE. 

214 KSX = P<6» 

216 PWR10X = 10.**(KSX-6) 

218 FX= P( 7>*PWR10X 
220 F=FX 
C 

222 1F(N0D( lFIX(P(5))t2).E0.1) GO TO 232 
224 TONLY=.TRUE. 

226 DX =P(8l*PWR10X 

232 IF(KSWI .E0.21N=NPTS 
11233 ILIM=N 

233 CALL PISTUGIXI 
IFCOX.EO.C.) GO TO 700 

234 FX=F 

248 IFLAB3=IFLAB3-KF0 
C 

250 IF(KSW64.E0.0)G0 TO 264 
252 KOUTX=-KPWR 
256 F10X=10.**KPWR 
260 WRITE (6,502) KOUTX 
C 

264 DO 278 1=1,11 

266 TEMP = FY+FLOATC I-1)*DY*10. 

268 ATEMP= ABS(TEMP) 

270 IF (ATEMP.LT.l.E-7) TEMP = 0. 

272 IF I ATEMP.GE.1.E+7)LAB0UT=2 
278 YLABEL( I )=TEMP 
300 KSYLAB =1 

302 WRITE (6,FYLAB) ( YLABE L ( I ) ,I =1 ,1 1 ) 
304 GO TO ( 306,700) , KSYLAB 
C 

306 KSYLAB =2 
310 LCTR=0 
NCTR=1 
K0UT*1 
KOUIT= 1 
C 

320 IF(MOO(LCTR,M) )328,322,328 
322 AFILL= XGRIO 
324 XGL =.TRUE. 

GO TO 330 
328 XGL =. FALSE. 

AF1LL=BLANK 
330 DO 332 1=2,104 
332 A( I I = AFILL 
334 DO 336 1=2, 104, NY 

336 A( I )= YGRID 
A(l) =BLANK 

338 GO TO ( 340,400),KOUT 
C 

340 XMIN=1.E15 
11340 IMIN=1 
11342 DO 11350 1=1,ILIM 
11344 IF( XMIN-X( 11)11350,11350,11346 
11346 XMIN=X( 1 ) 

11348 IMIN=I 
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11350 CONTINUE 

341 KX =iXf ININ)-FX)/OX +.5 

342 IF(KX-LCTR»630, 350*600 
350 LS=.TRUE« 

11352 X( IMIN 1-1.E15 
372 J*1 

374 IFIKSWS.EO.O) 60 TO 380 
376 ELS( J )-PI IMIN + 11) 

378 J= J + 1 

380 IFIKSW64.E0.0) GO TO 11400 
00382 ELS< J I=XMIN/F10X 
11400 DO 370 IM*1,KTL 
11402 LL = IMIN+I IM -1)*NPTS 
11404 KY»( Y(LL 1-FY1/DY+.5 
11420 IFIKSWI.E0.2) GO TO 11440 
11422 IK=0 

KLAST»2*KN+1 

11424 DO 11430 IL*3*KLAST*2 
11426 IK = IK+K( ILI 

11428 IFfIK-IMIN) 11430.11436,11436 

11430 CONTINUE 

11432 LA80UT*5 

11434 GO TO 700 

11436 KPC*K( IL+ll 

11438 GO TO 353 

11440 KPC=K(2»IM+2) 

353 KYL » KY+2 

354 IFIKY.lt. 0) 60 TO 360 
356 IFIKV.GT.IOIIGO TO 364 
358 GO TO 370 

360 KYL«1 
362 GO TO 366 
364 KYL=104 
366 TPC*RMARK 
370 A(KYL) *TPC 

386 IF(NCTR.GE.ILIM) GO TO 392 
388 NCTR*NCTR*1 
390 GO TO 340 
C 

392 KOUT = 2 
394 M* 10 
C 

400 IFIXGL.AND. fMOOlLCTR.lOl.EO.O) lKOUIT-2 

C 

600 WRITE f6*504)(A(I 1,1-1.1041 
602 IFIK864.NE.0) GO TO 620 
C 

604 1F(H0D(LCTR, 10)1614, 606, 614 
606 XLABEL *FX+FLOATILCTR)*DX 
608 TEMP sABSIXLABEL) 

610 IFITEMP.6E.1.E+7)LAB0UT=2 

612 IF! TEMP.LT.1.E-7)XLABEL=0. 

613 WRITE! 6, FLAB) XLABEL 
C 

614 LCTR = LCTR+1 

616 GO TO ( 320«302),KOUIT 
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c 


620 IF( .NOT. LSI GO TO 604 
62? LS = .FALSE. 

623 KT=K 864-1 

624 WRITE (6,FLSi(ELSf ll.lslvKTl 
626 GO TO 604 

630 L ABOUT 4 

700 GO TO (710* 702. 706,704.712) .LABOUT 
702 WRITE (6,506) 

GO TO 720 

704 WRITE) 6. 520) NCTR 

706 WRITE(6, 508)(X( 1),Y(I),I»1,2),K(1) ,K(2).K(3) 
708 GO TO 720 
712 WRITE(6,512) 

710 WRITE) 6, 510) 

720 RETURN 

500 FORMAT) 2HPT) 

502 FORMAT) 2HPT,13X,3HX*E, 12) 

504 FORMAT) 2HP ,26X,104A1) 

506 FORMAT) 2HPL.3X.10HBA0 LABELS ) 

508 FORMAT) 2HPL.5H N.G. .4G20. 8,31 7 ) 

510 FORMAT) 2HPL) 

512 FORMAT) 2HPL, 16HERR0R IN K ARRAY ) 

520 FORMAT) 2HPL.18HX OUT OF ORDER. Ia,I5) 

END 
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D. PLOTMY Dictionary 


AFILL 

ATEMP 

BLANK 

DX 

DY 

ELS 

F 

FlOX 

FLAB 


Array in which each line of the plot proper is constructed. It is ini- 
tialized to 104 blanks or 104 minuses (the horizontal grid-line charac- 
ter). The vertical grid-line characters (I's) are then placed in posi- 
tion 2 and every NY^^ position thereafter. Finally, the plotting char- 
acters (if any) are positioned. 

Temporary storage for the character (1 or -) with which the A array 
is initialized. 

Temporary storage for the absolute value of each y grid-label . 

BCD representation for a blank printing character. It is used to ini- 
tialize the A array for lines that are not x grid-lines. 

Vertical scale factor (per printing line). It is either calculated by 
PISTUG based on the range and length of the X array or it is calcu- 
lated from values supplied by the user of option 16 in P(6) and P(8) as 
DX = P(8)- It must be positive. Also see DX, COMMON 
/JOLO/ - Dictionary, page 53. 

Horizontal scale factor (per print position). It is either calculated by 
PISTUG based on the range and length of the Y array, or it is calcu- 
lated from values supplied by the user of option 32 in P(9) and P(ll) as 
DY = P(ll)*10^^^^‘®. Also see DX, COMMON/JOLO - Dictionary, 
page 53. 

Array in which the information (the first 18 characters of each line) to 
be printed at the left of the plot is assembled. See FLS on the follow- 
ing page. 

See COMMON/JOLO/ - Dictionary, page 53. 

IqKPWRX plox is the divisor of all values in the X array before 
they are printed on the left of the plot when option 64 is used. 

Variable FORMAT array with which each x grid-label is written. It 
is of the form 


(2HP+, 18X, F9.d) 

where 0 < d < 6. A value of d is inserted by the program (as KFD) 
after the call to PISTUG to process the X array. See KFD, page 53. 
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FLAB3 Initial value for FLAB(3). 


Variable FORMAT array with which the output to the left of the x 
grid-labels is written out if option 8 or option 64 is used. It is con- 
structed during execution and becomes one of three different state- 


ments. 


If 8 only (2HP+, 12X, A6) 

If 64 only (2HP+, 12X, F6.3) 
If 8 and 64 (2HP+, A12, F6.3) 


Used to write: 

One word of P array 

The modified x coordinate 

One word of P array and the 
modified x coordinate 


FLS264 Used to fill FLS(2) if option 64. 

FLS28 Used to fill FLS(2) if option 8. 

FLS2B Used to fill FLS(2) if both option 8 and option 64 are used. 

FLS38 Used to fill FLS(3) if option 8. 

FLS3B Used to fill FLS(3) if option 64 or both 8 and 64 are used. 

FORY See COMMON/ JO LO/ - Dictionary, page 53. 

FX Starting value for the vertical scale. It is either calculated by PISTUG 

based on the range, length, and XI (see p. 52) of the X array, or it is 
calculated from values supplied by the user of option 16 in P(6) and P(7) 
as FX = P(7)-10^^®^‘®. 

FY Starting value for the horizontal scale. It is either calculated by 

PISTUG based on the range, length, and XI (see p. 52) of the Y array, 
or it is calculated from values supplied by the user of option 32 in P(9) 
and P(10) as FY = P(10)- 

FYLAB Variable FORMAT array with which the line of y grid-labels is 
written. It is of the form 

(2HP , 20X, llFlO.d) 

where 0 < d < 6. A value of d is inserted by the program (as KFD) 
after the call to PISTUG to process the Y array. See KFD, page 53. 


FYLAB3 Initial value for FYLAB(3). 
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I 

IFLAB3 

II 

IK 

IL 

lUM 

IM 

IMIN 

lYLAB 

J 

K 

K864 

KFD 

KLAST 

KN 

KOBE 


Index of many loops. 

Equivalent to FLAB(3). 

Inner index of double DO loop that initializes the Y and K arrays for 
DUPY. 

Counter used as a subscript in the loop that associates the correct plot- 
ting character with the current value of IMIN for DUPY and NO DUP. 

Index of the loop that finds the correct plotting character for DUPY and 
NO DUP. 

Number of executions of the loop that finds XMIN and IMIN. 

Index of the loop that calculates KY, finds the associated plotting char- 
acter, and positions it in the A array. 

Value of I for the X(I) that is the smallest value in the X array at any 
particular time. 

Equivalent to FYLAB(3) . 

Subscript used when filling ELS array. 

Array name, the third argument in the call list of this subroutine. The 
first element is the sum of the option numbers being used; the second is 
the number of curves. The remaining elements of the array are pre- 
scribed by the Variation and the options being used. 

Sum of KSW8 and KSW64. It has the possible values 0, 1, 2, 3. 

See COMMON/JO LO/ - Dictionary, page 53. 

Limit of the loop that finds the correct plotting character for DUPY and 
NO DUP. 

Number of curves to be plotted. It is supplied by the programmer in 
K(2). 

Sum of the option numbers. It is supplied by the programmer in K(l). 
Branching on single bit positions using MASKl , MASK2, MASK4, 
MASK8, MASK16, MASK32, and MASK64 serves to identify the options 
being used. 
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KOUT 


Switch initialized to 1 and set to 2 after all points have been processed. 


KOUTX Value in the special heading of option 64 for the x coordinates. It is the 
negative of KPWR of the X array. 

KPC Equivalent to TPC , Temporary storage for a plotting character before 
it is positioned in the A array. 

KPCSTD Array containing the six standard plotting characters (* + 0 X = o). 

KPWR See COMMON/JO LO/ - Dictionary, page 53. 

KQUIT Switch initialized to 1. It is set to 2 just before the last line (a labeled 
X grid-line) of the plot is printed. 

KSW8 Switch initialized to zero and set to 1 if option 8 is used. 

KSW64 See COMMON/JOLO/ - Dictionary, page 53. 

KSWI Switch initialized to 1 and set to 2 if Variation I (DUPX) is used. 

KSWII Switch initialized to 1 and set to 2 if Variation II (DUPY) is used. 

KSX Fixed point form of the value supplied in P(6) by the programmer using 

option 16. P(7) and P(8) are multiplied by to get FX and DX. 

KSY Fixed point form of the value supplied in P(9) by the programmer using 

option 32. P(10) and P(ll) are multiplied by to get FY and 

DY. 

KSYLAB Switch initialized to 1 and set to 2 after the first line of y grid-labels 
has been written. 

KT Value of K864 - 1. It is the number of words written to the left of the 

plot; 1 if option 8 or option 64, 2 if both option 8 and option 64. 

KTIMES Value of KN - 1 . The number of executions of the loop that stuffs the 
Y and K arrays for DUPY. 

KTL Number of executions of the loop whose index is IM. It is 1 for DUPY 
and NO DUP and KN for DUPX. 

KX ’’Line number” of a particular point. It is computed from 

|(x. - FX)/D3^ rounded. The first x grid-line is ’’line number” zero. 
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KY 

KYL 

L 

LABOUT 

LCTR 

LL 

LS 

M 

MASKl 

MASK2 

MASK4 

MASKS 

MASK16 

MASK32 

MASK64 

MM 


"Printing position" of a particular point. It is computed from 

j(y. - FY)/DY rounded. The first y grid-line is in "printing position" 

zero - actually the second element of the A array. 

Value of KY + 2 if the y-coordinate falls on the plot; -1, if KY is nega- 
tive; 104, if KY> 101. 

Subscript computed and used in the loop that stuffs the Y and K arrays 
for DUPY. 

See COMMON/JOLO/ - Dictionary, page 53. 

Line counter. It is set to zero for the first x grid-line and is stepped by 
1 after each line is printed. 

Subscript computed and used in the loop that computes each KY and 
finds its associated plotting character. 

Logical variable set to TRUE whenever a point has been found for the 
current line. If option 8 or 64 is in effect, LS is used to control a 
branch (to print on the left side of the plot) and is set to FALSE imme- 
diately after the branch. 

Frequency of the x grid-lines. It is either set equal to 10 or supplied in 
P(3) by the programmer if option 2 is being used. If the programmer 
specifies zero, M is set to 1000, and only the first and last x grid-lines 
will be printed. 


>See KODE, page 41. 


Subscript computed and used in the loop that stuffs the Y and K arrays 
if DUPY is used. 
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N 

NCTR 

NPTS 

NPTST 

NY 

P 

PWRIOX 

PWRIOY 

RMARK 

STUG 

TEMP 

TUNX 

TONLY 

TPC 

X 

XGL 

XGRID 


See COMMON/JOLO/ - Dictionary, page 54. 

Counter of the number of points that have been processed. It is ini- 
tialized to 1 and tested for equality to ILIM. It is increased by 1 each 
time the test fails. 

Number of points in the X array (or Y array) supplied by the program- 
mer using DUPX (or DUPY) in K(3), 

Sum of the numbers of points in all the curves (supplied by the pro- 
grammer in K(3), K(5), . . . K(2 * KN + 1)) for NO DUP. It is KN 
multiplied by NPTS for DUPY. 

Frequency of the y grid-lines. It is either set to 10 or supplied in P(4) 
by the programmer using option 4. If it is supplied as zero, it is re- 
placed by 1000 and only the left-most y grid-line will be printed. 

Array name, the fourth argument in the call list of this subroutine. Its 
contents are prescribed by the Variation and the options being used. 

Value of 10^®^"®. 

Value of io^®Y-6 

Record mark character (+) . 

See COMMON/JOLO/ - Dictionary, page 54. 

Temporary storage. 

See COMMON/JOLO/ - Dictionary, page 55. 

See COMMON/JOLO/ - Dictionary, page 55. 

Current plotting character. Equivalent to KPC. 

Array name, the first argument in the call list of this subroutine, con- 
taining the values of the variable to be plotted on the vertical scale 
(down the page). The contents are destroyed during execution. 

Logical variable set to TRUE whenever the value of LCTR is evenly 
divisible by M. 

X grid-line character - a minus sign. 
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XLABEL 

XMIN 

XYX 

Y 

YGRID 

YLABEL 


Current value of the x label. It is computed only for every tenth line 
from XLABEL = FX + LCTR DX. 

Smallest value in the X array at any particular time. 

See COMMON/JOIXD/ - Dictionary, page 55. 

Array name, the second argument in the call list of this subroutine, 
containing the values of the variable to be plotted on the horizontal 
scale (across the page) . 

y grid-line character, the digit 1. 

Array holding the 11 labels for the y axis. 
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m. PISTUG 


A. Description 

This routine is internal to the plotting system and is used by PIXDTXY and 
PLOTMY. The calling statement is CALL PISTUG (ARRAY). ARRAY is an array in 
any order. All other information is transmitted through 

COMMON/JOLO/N, F, DX, XYX, FORY, STUG, LABOUT, TONLY, KSW64, KPWR, KFD, TUNX 

PISTUG uses the minimum and maximum values of the array to compute the total 
range. When scaling parameters are being computed for the y scale, 101 print posi- 
tions are available to cover the range. For the x scale, an arbitrarily chosen number 
of lines is assigned. In either case, the scale factor computed is always a value of 
DxlO^, where D is a member of the set [2, 2.5, 5, lO]. 

PISTUG uses a function subprogram KHAR(V) which finds the integer that is the 
characteristic of the logarithm to the base 10 of the floating point variable V. 
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B. Program Listing 


SUBROUTINE P I STUG ( ARRA Y» 

DIMENSION ARRAY! 1 ) 

COMMON /J0L0/N,F,DX.XYX,F0RY,STUG,LAB0UT»T0NLY,KSM64,KPWR,KFD.TLINX 
LOGICAL XYXtFORY, STUG«TONLY 
126 XI s ARRAY! 1) 

128 IF(XYX) GO TO 133 
130 DO 132 J * 2.N 

132 XI - AMINK XltARRAYUn 

133 IF(STUG) X1=F 

134 XN = 0.0 

136 DO 146 J * l.N 

138 OIF 3 ABS(Xl-ARRAYIJ) ) 

140 IF! OIF.LE.XNIGO TO 146 
142 XN*DIF 
144 I HOLD* J 

146 CONTINUE 

147 XN = ARRAY! IHOLO) 

148 IFIKSW64.E0.2) KPWR « KHAR ! AMA XI lABS ( XI ) .ABS ! XN) H 

149 IFITONLY) GO TO 240 

150 TLIN*101. 

152 IF! .NOT.FORY) TLIN * TLINX 
154 C5 * !XN-X1)/TLIN 
156 C6 » ABS!C5) 

158 IF!C6.E0.0.) GO TO 300 

159 K7 a KHAR(C6) 

160 C8 a 10. ♦♦K7 
162 C9 a C6/C8 

164 IF! I 2.5-C9).LE.O.O) GO TO 172 
166 0a2. 

168 IF! 1 2.0-C9).LE.0.) Da2.5 
170 GO TO 176 
172 Da5. 

174 IF! !5.-C9).LE.0.0) 0*10. 

176 Cll a d»C8 

178 DX a SIGN! Cll, C5) 

179 HUND a 100. ♦OX 

240 K7 a KHAR!ABS(OXn 
250 KFD a 0 

252 IF1K7) 260,270,254 
254 IF!K7.GE.5) LAB0UT=2 
256 GO TO 270 
260 KFO a 6 

262 IF!K7.LT.I-7) ) LABOUT = 2 
264 IF!K7.6T.!-6) ) KFO a -k7 
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270 IF(STUG) GO TO 230 
182 KC.12 = INTtABS«Xl)/Cin 
184 JJ = 1 

186 IF( XI I 188, 192,190 
188 JJ = 3 

190 IFCDX.LT.O.) JJ = JJ+1 

192 GO TO ( 193, 194, 195,196) ,JJ 

193 KC14 = KC12 

GO TO 198 

194 KC14 = KC12+1 

GO TO 198 

195 KC14 = -KC12-1 

GO TO 198 

196 KC14 = -KC12 

198 KC13 = M00IKC12,10) 

KC15 = KC12-KC13 

199 KC18 = KC15 

200 GO TO C 212, 202, 202, 210), JJ 
202 KC18 = KC18+10 

204 IF(KC13.NF.9) GO TO 210 
206 KC18 = KG 14 
20 8 GO TO 212 

210 IF1X1.lt. 0.)KC18 = -KC18 
212 F=C11*FL0AT(KG18) 

214 IFI .NOT.FORY) GO TO 230 
220 TEMP = F+HUND 
222 GO TO I 224, 228, 224, 228), JJ 
224 IFI TEMP. GE.XN) GO TO 230 
226 GO TO 229 

228 IF! TEMP.LE.XN) GO TO 230 

229 F=C114FL0AT(KC14> 

230 CONTINUE 

RETURN 
300 DX=0. 

I A80UT=3 
GO TO 230 
END 


FUNCTION KHAR(XMAX) 

KHAR = INTIAL0GlXMAX)/2. 302585+40. 0-40 

RETURN 

END 
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D. PISTUG Dictionary 

C5 Value of the maximum possible scale-factor. It is found by dividing 

the range XN - XI by TUN. 

C6 Absolute value of C5. It may not be zero. If it is, the error switch 

LABOUT is set to 3 and the plot is terminated. It is used to deter- 
mine C9. 

C8 Intermediate storage used in the calculation of C9 and DX. 

C9 When the maximum possible scale -factor C5 is expressed as x.xxxx 

multiplied by 10*^, C9 is the absolute value of the x.xxxx factor. It is 
used to select D, the next larger value from the set of permissible 
values [2, 2.5, 5, 10]. 

Cll Intermediate storage used in the calculation of DX and of F. 

D Smallest member of the set [2, 2.5, 5, lO] that is larger than C9. It 

is used when PISTUG is calculating the scaling parameter DX. 

DIF Absolute value of the distance between XI and the element furthest 

from XI . 

DX See COMMON/JOIXD/ - Dictionary, page 53. 

F See COMMON/JOLO/ - Dictionary, page 53. 

FORY See COMMON/JOLO/ - Dictionary, page 53. 

HUND One hundred times the scaling parameter DX. 

IHOLD Index of the element in the array which is furthest from XI. 

J Index of the search loop in which DIF and IHOLD are found. 

JJ Switch used to control the calculation of the starting value, F, 

according to the algebraic signs of XI and DX. 

K7 Characteristic of the logj^Q of the absolute value of the scale-factor 

DX. It controls the value of KFD, the number of decimal places in 
the grid-label, as follows: 
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KC12 

KC13 

KC14 

KC15 

KC18 

KFD 

KHAR 

KPWR 

KSW64 

LABOUT 

N 

STUG 

TEMP 

TUN 

TUNX 

TONLY 


KFD = 0 if K7 is + or zero. 

KFD = -K7 if K7 is negative and greater than -6. 

KFD = +6 if K7 is negative and equal to or less than -6. 

If K7 is greater than 4 or less than -7, the error switch LABOUT is 
set to 2 to cause an error message to be printed below the plot. 

The symbol K7 was inadvertently also used for intermediate storage of 
the characteristic of logj^Q C5 during the calculation of C9. 


I Intermediate storage used in the calculation of the "best possible" 
I value of F. 


See COMMON/JOLO/ - Dictionary, page 53. 

Function subprogram (see p. 46). 

See COMMON/JOLO/ - Dictionary, page 53. 

See COMMON/JOLO/ - Dictionary, page 53. 

See COMMON/JOLO/ - Dictionary, page 53. 

See COMMON/JOLO/ - Dictionary, page 54. 

See COMMON/JOLO/ - Dictionary, page 54. 

Temporary storage. 

When PISTUG is scaling a Y array, TUN is set equal to 101. When 
an X array is being scaled, TUN is supplied by the calling routine as 
TUNX. TUN is used to compute C5. 

See COMMON/JOLO/ - Dictionary, page 55. 

See COMMON/JOLO/ - Dictionary, page 55. 


PISTUG 


51 


Minimum or maximum value of the array being scaled. It is set equal 
to the first element of the array when PISTUG is scaling an X array 
for P1X)TXY, otherwise it contains the minimum value of the array 
after a search. If STUG is TRUE, XI is set equal to the scaling- 
parameter F supplied by the calling routine. XI is used to compute 
C5 and KPWR, and several times in the calculation of F. 

Minimum or maximum value of the array being scaled. It is set equal 
to zero and replaced by the element furthest away from XI . It is used 
to compute C5 and in the calculation of KPWR. 

See COMMON/JOLO/ - Dictionary, page 55. 
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DX 

F 

FORY 

KFD 

KPWR 

KSW64 

LABOUT 

1 


Scaling parameter for the scale factor for either scale. It may be cal- 
culated by PISTUG or supplied indirectly by the user. For PLOTMY, 
DX for the x scale must be positive. 

Scaling parameter for the starting -value for either scale. It may be 
calculated by PISTUG or supplied indirectly by the user. 

Logical switch set to TRUE by PIGTXY or PIX)TMY only if the call 
to PISTUG is for the y scale . 

Number of decimal places required in the grid-label for the array 
being processed by PISTUG. It is used in PLOTXY and PLOTMY to 
modify the variable FORMATS (FLAB, FYLAB) used to write the 
grid-labels. This modification is made by subtracting KFD from a 
word in the FORMAT array containing [blank, blank, F9.0] . Sub- 
traction rather than addition is necessary because the representation 
of a blank character in the 7094 is 110000, which is negative. 

Characteristic of the logjQ of the maximum of the absolute values of 
XI and XN. (The minimum and maximum of the array being pro- 
cessed by PISTUG. ) It is calculated by PISTUG only if option 64 is 
being used. It is used by PIOTXY and PLOTMY to prepare the coor- 
dinates being printed at the left of the plot so they will fit a F6 . 3 
FORMAT specification. A special heading is written displaying the 
conversion factors. 

Variable initialized to zero by PLOTXY and PLOTMY and set to 2 if 
option 64 is used. It is used in PISTUG to control calculation of 
KPWR and in PLOTXY and PIOTMY to control printing of coordinates 
at the left of the plot. 

This is a computed-go-to index used to control branching to normal 
termination or to various error terminations. 

Normal ending of a plot; a blank PL record. 
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Error message BAD LABELS is written below the plot. This is an 
indication that on e or more x or y grid-labels are either too large 
or too small to be printed correctly in the assigned label field 
(F9. d - where 0 < d < 6) . This is checked by PISTUG for F and by 
PLOTXY and PLOTMY for successive labels. The relative positions 
of the plotted points are always correct; the labels may or may not be. 

The plot is terminated immediately if: 

(a) The range of an array being scaled by PISTUG is found to be zero. 
The error output is 

PLOTXY N. G. followed by X(l), Y(l), X(2), Y(2), K(l), P(l) 
PLOTMY N. G. followed by X(l), Y(l), X(2), Y(2), K(l), K(2), K(3) 

(b) The values in the X array are not in monotonic order when 
PLOTXY is used. The error output is X OUT OF ORDER followed 
by the message in (a) . 

The plot is terminated immediately if option 16 is used in PLOTXY or 
PLOTMY and a value of x is found outside the user-specified starting- 
value (Rewind Carriage Error) . The error output is the same as in (b) . 

If the search in PLOTMY for the plotting character fails, the plot is 
terminated and the error message ERROR IN K ARRAY is written. 

Number of points in the array sent to PISTUG to be scaled. If the call 
is from PLOTXY, N has been supplied by the user in P(l). If the call 
is from PLOTMY: 

Variation I - N is the number of points in the X array (supplied 
by the user in K(3)) . 

Variation II - N is the number of points in the X array. It is cal- 
culated by multiplying the number of points in the Y array 
(supplied by the user in K(3)) by the number of curves, KN 
(supplied by the user in K(2)). 

Variation HI - N is the total number of points to be plotted and is 
calculated by adding the number of points for each curve 
(supplied by the user in K(3), K(5), K(7), etc). 



STUG 


TUNX 

TONLY 

XYX 


Logical switch set to TRUE by PLOTXY or PLOTMY when PISTUG is 
not required to calculate F. 


If PISTUG is scaling an X array, the calling routine must supply the 
number of lines into which the plot must be fitted. The empirically 


chosen function of N, 


TUNX = 55^— + ij 


35 points onto one page. 


fits most plots of less than 


Logical switch set to TRUE by PLOTXY or PLOTMY if PISTUG is not 
required to calculate F or D. 

Logical switch set to TRUE only if the call to PISTUG is for the x 
scale from PLOTXY. 
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IV. SORTXY 


A. Description 


CALL SORTXY (V, W, NPTS) 

This subroutine rearranges the NPTS elements of the V array in order of in- 
creasing size. The elements of the W array are moved to maintain the original pair- 
relation; that is, if the fifth element of the V array is moved to the first position of 
V, the fifth element of W is moved to the first position of W. 

B. Program Listing 


SUBROUTINE SOR T XY( X ♦ Y ,NP TS) 
DIMENSION X( I ) ♦ Y( 1 » 

100 N=NPTS 
in? NN=N-1 
104 00 '40 KT=1.NN 
XMIn 
JAD=K . 

JKL=KT rl 

11? 00 120 JK=JKL.N 

114 IF (XMIN-XUK)J 120,120,116 

116 XMIN=X(JK) 

118 JAD=JK 
120 CONTINUE 
12? YMIN=Y(JAD» 

X(JAD)= X(KT» 

Y(JA0)= Y(KT) 

X1KT)= XMIN 
Y(KT)= YMIN 
140 r.ONTINU 
RETURN 
END 
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V. SCALE 


A. Description 


CALL SCALE (NPTS, A, KRSTR) 

This subroutine finds the largest absolute value of the NPTS elements of A, and 
computes the characteristic of its logjQ. If the characteristic K is -2 < K < 4, 
KRSTR is set to zero and control returns to the calling program. If K > 4 or <-2, 
each element of A is multiplied by a power (KRSTR) of 10 to transform the array to 
suit PLOTXY and PLOTMY. 

KRSTR is returned to enable the user to ReSToRe the array or record how it has 
been altered. 

B. Program Listing 


SUBRUUTINJE SCALECNPTS, X,KRSTRI 
DIME.NSIOM X(1) 

100 N=NPTS 
102 XMAX=X( 1) 

104 DO 1U6 J = l,'4 

106 XMAX=AMAX1(ABS(XMAX»,ABS(X(J)) J 

108 KHAR = liMT(AL0G(XMAX)/2.302585+40. »-40 

116 IF( (4-KHARJ<'(KHAR+2»)120,118,118 

118 KPWR=0 

119 GO TO 130 

120 KPWR=3-KHAR 
122 FACr=lC.4*KPWR 
126 00 128 J=1,N 
128 X( J )=X(J )#FAtT 
130 KRSTR=KPWR 

140 RETURN 
END 
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VI. INPUT-OUTPUT CONSIDERATIONS 


A. Discussion 

The printer plotting system used at Lewis operates under a supervisory system 
which controls all input-output. All BCD records (limited to 132 characters) written 
with a WRITE (6,###) statement are transferred from the 7094 core to the 7044 core 
to the disk and later from the disk to the 7044 core to one of two 1403 printers. 

The 7044 program which controls printing must handle records written by the 
plotting routines somewhat differently than ordinary output records. For this reason 
the first character of every record written by PLOTXY and PLOTMY is a P. The 
P is recognized by the supervisory program and automatic page-spacing and normal 
first-record-character printing controls are disabled. Further processing is based 
on the second character in each such P record. The flow chart (p. 60) displays the 
necessary logic. 

For your convenience in making input-output changes the WRITE and FORMAT 
statement numbers for PLOTXY and PLOTMY are listed below. 

WRITE FORMAT 

100 500 Writes an otherwise blank PT record. May 

be the only one and will cause a skip to new 
page, or it will be the last one and will cause 
a blank line. 

260 502 Writes a special heading to the left of the 

plot if option 64 is used. 

302 FYLAB This FORMAT statement is used to write the 

lines of y grid- labels and is of the general 
form (2HP , 20X, llFlO.d). It is an array to 
permit setting of the value of "d" (in the F 
specification) to suit the range of the Y array. 
See KFD, page 53. 

600 504 Writes each line of plot after grid-line char- 

acters and plotting characters have been po- 
sitioned. 


I 


58 


WRITE 


FORMAT 


624 


FLS 


PLOTXY 

PLOTMY— 

613 FLAB 




702 

506 


704 

520 

► 

706 

508 

J 

710 

510 



PLOTMY ONLY 
712 512 


This FORMAT statement is used to print at 
the left of the plot. It is constructed during 
execution and becomes one of three different 
statements ; 

If option 8 (2HP+, 12X, A6). 

r If option 64 (2HP+, 6X, 2F6.3). 

I^If options 8 and 64 (2HP+, A6, 2F6.3). 

pf option 64 (2HP+, 12X, F6.3). 

I^If options 8 and 64 (2HP+, A 12, F6.3), 

This FORMAT statement is used to write 
each X grid-label and is of the general form 
(2HP+, 18X, F9.d). It is an array to permit 
setting the value of *'d" (in the F specifica- 
tion) to suit the range of the x array. See 
KFD, page 53. 


These are error messages. 


This FORMAT statement is used to write an 
otherwise blank PL record. It may be the 
only such record and one is necessary to re- 
set switches in the printing routine. If there 
are more than one, this will be the first and 
it will cause a blank line below the lower line 
of y grid-labels. 


This is an error message if the search for 
the associated plotting character fails. See 
page 14. 
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B. Logical Flow Chart. 


Any output 
record is 
ready 


Normal print 
procedures 


character 


3witch> 


Set® 
to ON 


Skip to new 
page. Set ® 
to OFF. 

Set ★ to 
FIRST. 


characteiV 
V T? / 


switch^ 


FIRST 


Skip to new page. 

Set ★ to NOT-FIRST. 


NOT-FIRST 


character 
?/ 



Replace second char- 
acter with blank and 
print normally. 


Replace second char- 
acter with blank and 
print on same line as 
previous record (space 
suppress). 


Second 

\characten 


Set ★ to 
FIRST 


Switch ® must be loaded set to OFF, but need not be reset to OFF for each user nor for each plot. 
Switch ★ must be loaded set to FIRST, but need not be reset to FIRST for each user nor for each plot. 







SAMPLE PLOTS 



X»EXP(-3.6iJ-y*Y/4. 5) 
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FOUR VALUES QF THE OEPENUENT VARIABLE (ACROSS) 
FOR EACH VALUE OF THE INDEPENUENT VARIABLE (DOWN) 
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Figure 3. - This plot displays the use of PLOTMY - Variation I, KODE • 
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